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Abstract 

The probability of a random polygon (or a ring polymer) having a knot type K 
should depend on the complexity of the knot K. Through computer simulation using 
knot invariants, we show that the knotting probability decreases exponentially with 
respect to knot complexity. Here we assume that some aspects of knot complexity 
are expressed by the minimal crossing number C and the aspect ratio p of the tube 
length to the diameter of the ideal knot of K, which is a tubular representation of K 
in its maximally inflated state. 

1. Introduction 

Various species of knotted polymers have been synthesized and observed in chemistry and 
biology in the last two decades. Once a ring polymer is formed, its topological state is 
unique and invariant. The topological constraint on the ring polymer should be nontrivial. 
It may restrict the available degrees of freedom in the configuration space of the polymer, 
to a great extent. Consequently, it may lead to a large entropic reduction, which is related 
to the probability of random knotting, as wee shall see shortly. 
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For a knot we define the knotting probability Pk{N) by the probability that the 
topology of a random polygon with nodes is given by the knot K. If a ring polymer is 
under the topological constraint of the knot K, then the decrease of the polymer entropy is 
given by ASk = —ks log Pk{N), where Pk{N) is the ratio of the volume of the configuration 
space under the topological constraint to that of no topological constraint. The knotting 
probability Pk{N) should also correspond to the probability that a ring polymer of 
Kuhn units have the knot K when it is closed randomly during its synthesis. The knotting 
probabilities have been measured as the fractions of knotted species of circular DNAsii 

Let us now discuss how to express the complexity of knots. We could classify knots 
completely, if we might know all the topological properties that are invariant under any 
continuous deformation of the spatial configurations. The topology of a given polygon can 
be effectively detected by calculating some topological invariants such as the Alexander 
polynomial A/<(t) and the Vassiliev-type invariants Vn{K). Although the invariants are 
practically useful for computer simulationsi'i, it is not easy to derive any explicit topological 
properties or meanings from them. Let us consider the minimal number of crossing points 
in the knot diagram of a knot K. We donate it by C, or \K\ for the knot K. The minimal 
crossing number C should be a measure on the complexity of knots. The number C is useful 
in studying statistical or dynamical properties of knotted ring polymer There is some 
arguments on the mean-square radius of gyration of knotted ring polymers with respect to 
However, C is rather weak as a topological invariant. The number of knots that have 
the same number C increases rapidly: there are 166 primes knots which have 10 crossings. 

Recently the concept of ideal knots has attracted much interest. One of the most 
ideal (or elegant) geometric representations of a knot should be given by such a closed tube 
with uniform diameter that gives the largest ratio of the diameter to the tube length. We 
call such geometric representations ideal knots, briefly. For a ring polymer with a knot type 
K, Grosberg et.a/0 discussed a topological invariant p of a knot K, which is defined by the 
aspect ratio of the length (L) to the diameter (d) of such a ring polymer of knot K that 
is maximally inflated, i.e., an ideal knot. Here, p is given by L/d. We also denote it by 



Pk- The value pk should be a measure of knot complexity. It may be more powerful than 
C, since different knots should have different values of p, in general. Katritch et. al have 
obtained ideal knots for 42 different knotsi'0. There is a linear relation between the average 
crossing numbers of ideal knots and their p values.S0. The value p should be also useful for 
describing flexible DNA knots in thermal equilibrium!. Furthermore, the value p should be 
useful for statistical or dynamical researches of knotted ring polymers.liiS 

The A^-dependence of the knotting probability has been studied 
through simulations,i'&ii and it is found that the probability of the unknot (the trivial 
knot) decreases exponentially with respect to A^: 

Po(iV) = Coexp(-Ar/Ar,), (1) 

where Co and are fitting parameters. For some nontrivial knots (3i, 4i, 5i, 62), knot- 
ting probabilities have been evaluated numerically for several different models of random 
polygons and self-avoiding polygons. Through the simulations using the Vassiliev-type 
invariants, it is found that the probability Px{N) as a function of can be expressed as 

Here Ck, Nk and rriK are fitting parameters to be determined from the numerical results. 
The expressions (|l]) and (|^) should correspond to the asymptotic expansion of renormal- 
ization group arguments. Numerically we see that the estimates of should be given by 
almost the same value for any knot and therefore the N^^s are almost equal to A"c, 

which depends on the model. We also observe that the value mx of a knot K should be 
universal for the different models^. 

In this paper, we discuss how the knotting probability Pk{N) of a knot K should depend 
on its complexity while A^ being fixed, or in short, the knot dependence of the normalization 
constant Ck- Evaluating the knotting probabilities of several prime knots for Gaussian ran- 
dom polygons, we observe a rough tendency that the amplitude Ck decreases exponentially 
with respect to p. The numerical result seems to be favorable to Grosberg's conjecturelii 
that the probability Pk{N) as a function of p should be given by 



Pk{N) ~ exp(-iV/iV, - sp). (3) 

Here s is a constant. At this stage, however, we could not judge whether the conjecture 
be vahd, since the data points scatter outside the range of statistical errors. On the other 
hand, we show another explicit statistical behavior for a version of the knotting probability. 
Let us define the average knotting probability Pave{N, C) by 

PaUN,C)= Pk{N)/Ac. (4) 

K: \K\=C 

Here the sum is over such knots that have the same minimal crossing number C, and Ac 
denotes the number of prime knots which have the same minimal crossing number C . For 
instance, we have ^3 = Ai=l, A^ = 2 and Aq = 3. Then, we shall see from the data that 
the average knotting probability decreases exponentially with respect to C. Furthermore, if 
we consider the average of the p's over such knots that have the same C 

{p)= E Pk/Ac, (5) 

K:\K\=C 

then we see that the average knotting probability also decays exponentially with respect to 
the average (p). 

2. The method of simulations 

Using the conditional probabilitylli, we construct a large number of Gaussian random poly- 
gons, say M polygons, of nodes for = 300, 500 and 1000. Then, the knotting probability 
of a knot K is evaluated by Pk{N) = M^/M. Here Mx is the number of polygons with the 
knot K, and M is given by M = 10^ to each of the three numbers of A^. 

In order to detect the knot type of a give polygon, we employ three knot invariants: the 
determinant \Ax{t = — 1)| of a knot K, the Vassiliev-type invariants V2{K) and |f3(ir)| of 
the second and third degrees, respectively. We evaluate Mk, after enumerating the number 
of the polygons which have the same set of the values of the three invariants for the knot K. 
Using |f3(i^)|, we do not distinguish the chirality of the knot, i.e., the right-handed knots 



and the left-handed one^E. Furthermore, we do not consider six knots (89, 810, 811, 818, 
820 and 821) in any of the simulations in the paper. They have the same values of the three 
invariants as those of some composite knots. 

3. Results and Discussion 

The estimates of the average knotting probability Pave{N,C) are plotted in Fig. 1 against 
the minimal crossing number C, up to C = 8. It is clear in Fig. 1 that the average knotting 
probability Pave{N, C) decreases exponentially with respect to C. We remark that error bars 
correspond to one standard deviation in all the four Figures in the paper. 

Let us now discuss the knotting probability in terms of the p values. In Fig. 2, the 
knotting probabilities Pk{N) for some prime knots are plotted against px, where N is 
kept constant. We note that the p values of the 42 knots are listed in Table 1 of Ref.0, 
which are used in the paper. We see in Fig. 2 that there is a rough tendency that the 
knotting probability of a prime knot decrease exponentially with respect to the value p. 
The observation should be useful. However, it seems that there is no clear relation between 
the knotting probability Pk{N) and the aspect ratio p, since the data points of larger p 
values deviate from the possible regression line, considerably. Here we recall that error bars 
correspond to one standard deviation. 

Let us discuss the knotting probability for such knots that have the same minimal crossing 
number. For instance, there are two knots with 5 minimal crossings: 5i and 82. For Gaussian 
polygons, the knotting probability of 82 is always larger than that of 5i. This is consistent 
with the simulation of the cylindrical self-avoiding polygons^. Let us consider the three 
prime knots with C = 6. We observe that the knotting probabilities of 61 and 62 are almost 
the same, while that of 63 is always smaller than the other two. For prime knots with C = 7 
or 8, the data points are so close to each other that it is difficult to give any definite ranking 
on them. 

In terms of the average value {p), which is a function of C, the estimates of the average 
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knotting probability PaveiN,C) are expressed in Fig. 3. We clearly see the exponential 
decay of the average knotting probability Pave{N, C) with respect to (p). It is similar to Fig. 
1. This result shows that the entropy of a ring polymer with knot K decreases with respect 
to knot complexity expressed in terms of {p) or C. Here it is also suggested that {p) should 
be approximately linear to C. 

Let us discuss the iV-dependence of the knotting probability in terms of knot complexity. 
The ratio Pk{N)/ P^^^N) of a knot K against the value p is plotted in Fig. 4 for the three 
numbers of A^: N = 300, 500, 1000. Here we note that the trefoil knot (3i) is dominant 
among the nontrivial prime knots for the three iV's. We find again the rough tendency that 
the ratio Pk{N) / P^-^lN) decays exponentially with respect to knot complexity p. Moreover, 
for any knot K, the ratio Pk{N)/ P^^^lN) is given by almost the same value for the three 
numbers of A^, with respect to error bars as seen in Fig. 4. Thus, the ratios Pk{N) / P^^lN) 
are independent of A^. 

The above observation in Fig. 4 can be explained by using the fitting formula (|^). Let us 
assume that for a prime knot K, the exponent mx of eq. (0) should be given by almost the 
same value. Then, we have Pk{N)/ P^^^i^N) ~ Ck/C'^^, which is clearly independent of A^. 
Thus, in terms of the formula (^, the rough exponential decay of the knotting probability 
with respect to p is closely related to the knot complexity- dependence of the amplitude Ck- 
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FIGURES 

Fig. 1. Average knotting probability Pave{N,C) versus the minimal crossings C for 29 prime 
knots with upto C = 8. The line is given by Pavei^-, C) = Pavei^, 0) exp(— aC) with a = 1.16. 

Fig. 2. Knotting probability Pk{N) with N = 500 versus the aspect ratio p for 29 prime knots 
with upto C = 8. 

Fig. 3. Average knotting probability Pave{N,C) versus the average (p) for 29 prime knots with 
upto C = 8. The line is given by Pave{N, C) = Pave{N, 0) exp(-/? (p)) with /3 = 0.30. 

Fig. 4. The ratio of knotting probabilities Pk{N)/P3j^{N) versus the aspect ratio p for 
= 300, 500, 1000. 
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